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Abstract 

Anisotropy data analysis leaves a significant degeneracy between primeval spectral 
index (n s ) and cosmic opacity to CMB photons (r). Low-/ polarization measures, in 
principle, can remove it. We perform a likelihood analysis to see how cosmic variance 
possibly affects such a problem. We find that, for a sufficiently low noise level {c^ ix ) 
and if r is not negligibly low, the degeneracy is greatly reduced, while the residual 
impact of cosmic variance on n s and r determinations is under control. On the 
contrary, if a^ ix is too high, cosmic variance effects appear to be magnified. We apply 
general results to specific experiments and find that, if favorable conditions occur, 
it is possible that a 2— a detection of a lower limit on r is provided by the SPOrt 
experiment. Furthermore, if the PLANCK experiment will measure polarization 
with the expected precision, the error on low-/ harmonics is adequate to determine r, 
without significant magnification of the cosmic variance. This however indicates that 
high sensitivity might be more important than high resolution in r determinations. 
We also outline that a determination of r is critical to perform detailed analyses on 
the nature of dark energy and/or on the presence of primeval gravitational waves. 
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1 Introduction 



Data on cosmic microwave background (CMB) anisotropy and polarization 
can provide effective constraints on cosmological parameters. However, while 
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significant anisotropy observations are available (see, e.g., Smoot et al. 1992; 
de Bernardis et al. 2000; Hanany et al. 2000; Halverson et al. 2002; Sievers 
et al. 2002; Scott et al. 2002) and further anisotropy data are expected from 
experiments in progress, detailed data on the weaker polarization signals may 
become available only within a few years. In turn, dark energy (DE) param- 
eters or the amplitude of primeval gravitational waves (GW) are just mildly 
constrained by anisotropy data, while polarization would significantly depend 
on them. A peculiar situation then holds for the cosmic optical depth to CMB 
photons (r), due to reionization. Anisotropy data only mildly constrain it, 
while polarization data would provide more stringent constraints; of course, 
greater r values would ease its detection, but, in turn, greater r values would 
also ease the determination of DE and/or GW parameters. 

This paper tries to investigate which polarization experiment (s) can determine 
t. The success of an experiment, however, will also depend on the value of r 
itself. In particular, we shall consider the SPOrt experiment (Macculi et al. 
2000; Carretti et al. 2000; Peverini et al. 2001), planned to perform a (nearly) 
full-sky polarization measure aboard of the ISS in 2004, with an angular 
resolution similar to COBE. 

When only anisotropy data are considered, there is a degeneracy between 
the effects of varying r or the primeval spectral index n s (Jungman et al. 
1996; Zaldarriaga, Spergel & Seljak 1997; Eisenstein, Hu & Tegmark 1999; 
Efstathiou & Bond 1999). The most stringent limits on r, based on anisotropy 
data, are provided by Stompor et al. (2001), who used the outputs of recent 
balloon experiments to set r < 0.4 at the 2-a level, once n s is assumed 
< 1.2. (Such limits will be included in a number of figures here below and 
denominated Stompor limits). Full-sky anisotropy measures from space, like 
MAP 3 or PLANCK 4 , might improve the Stompor limits by a factor ~ 2. In 
turn, a residual uncertainty (~ 0.05) on the value of n s is related to our scarce 
knowledge of r. Of course, a more precise determination of n s would really be 
welcome, first of all to shed new light on the mechanism generating primeval 
fluctuations. Furthermore, within a generic inflationary model, only for n s < 1 
primeval GW may be expected to exist. 

The impact that a variation of r may have on CMB polarization can be seen 
from Fig. 1; here, for a spatially flat cosmological model (density parameters: 
Vt m = 0.35, fit = 0.05; Hubble constant in units of 100 km/s/Mpc: h = 0.65; 
n s — 1; dark energy is an ordinary cosmological constant) and for different 
values of r, we plot the dependence on / of the angular spectra Q 4 with 
A = T, E, B, X, standing for anisotropy, E- and S-mode polarization and 
anisotropy-polarization cross correlation respectively; only E-mode polariza- 
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tion is considered through this paper. Fig. 1 shows that the main dependence 
on r actually occurs at low I and for the polarization and cross-correlation 
spectra; comparatively, the r dependence of spectra, at I > 30, is milder. This 
point is important to devise the right polarization experiment, which should 
aim to great sensitivity, while, for r determination, high resolution data might 
not bear a direct relevance. Accordingly, the angular resolution of the SPOrt 
experiment does not prevent it from achieving cosmologically significant re- 
sults. 

Here, as well as in most recent literature, we assume that the whole cos- 
mic opacity arises from an almost complete ionization of hydrogen, after a 
suitable redshift. More complex reionization histories may have to be taken 
into account to discuss future data (Bruscoli, Ferrara & Scannapieco 2002; 
Venkatesan 2002), but would not modify general conclusions. 

For low-l spectral components, varying n s has a milder effect than varying 
r. If the whole angular spectrum is considered, other parameters, like fltot, 
Q m , Qb (total, matter and baryon density parameters), may be determined 
independently from low-/ features. Here we shall assume that they are already 
known from such high-/ anisotropy measures. 

Dealing with low Z's, a significant part of our analysis will be devoted to inspect 
cosmic variance. Its impact on C\ determination is known. Here, however, we 
are interested in measuring different quantities and a knowledge of C;'s is 
not even required. We shall explore this point by simulating and analyzing a 
large number of artificial data sets. Among other things, we shall see that the 
SPOrt experiment, although far from granting a stringent r determination, 
may improve the Stompor limit in a significant number of cases. 

In the next section we shall first debate which range of values can be expected 
for r. We shall then describe how artificial data sets are produced. In section 4 
we shall show the results of their analysis, in a number of different experimental 
conditions and considering different r values. In section 5, conclusions will be 
drawn. 



2 The cosmic opacity to CMB photons 

In the present cosmic epoch, although diffuse baryonic materials are almost 
completely ionized, the scattering time for CMB photons t s ~ 4.45-10 18 fi^ 1 /i~ 2 s 
(h: Hubble parameter in units of 100 km/s/Mpc) exceeds the Hubble time 
tu = 3.09 • 10 17 /i~ 1 s and the Universe is transparent to CMB photons. These 
figures can be extrapolated to the past, but, while t s oc a 3 (a: the scale factor), 
the dependence of t# on a depends on the cosmological model. In an expan- 
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sion regime dominated by non-relativistic matter, a fully ionized Universe is 
opaque to CMB at z > z op ~ 5 {fl b h)- 2 ^. Therefore, assuming Vl^h 2 < 0.022 
(standard BBNS limits; see, e.g., Dolgov 2002 for a recent review) and h ~ 0.5, 
we have z op > 64, safely above the expected reionization redshift z ri . No rea- 
sonable cosmological model allows a substantial reduction of such z op . 

On the contrary, the precise value of r(^C 1), does depend on the model and, 
in particular, on z ri . Recent data on high-z QSO's (Djorgovski et al. 2001; 
Becker et al. 2001) indicate that at least a fraction of neutral hydrogen is 
present, in the intergalactic medium, at z > 6. These very authors, however, 
warn us against an immediate conclusion that z„ ~ 6. Reionization can be 
a slow and/or patchy process and our present knowledge of primeval object 
formation can be consistent with an effective z ri ~ 6 or a significantly greater 
value. In Fig. 2 we show how r is related to z ri and model features, considering 
a wide set of cosmological models. For z ri < 15, values of r up to ~ 0.15-0.20 
are licit. 



3 The artificial data sets 

Artificial data sets will be built with reference to the features of the Sky Po- 
larization Observatory (SPOrt) experiment, selected by the European Space 
Agency (ESA) to fly on the International Space Station (ISS) in 2004, aiming 
to measure the linear polarization of the diffuse sky radiation at 22, 32 and 
90 GHz, using a set of polarimeters fed by corrugated horns with an angular 
resolution (FWHM) of 7°. At variance from other space experiments (MAP 
and PLANCK), SPOrt has been explicitly designed to measure the Stokes 
parameters Q and U and to minimize systematics and instrumental polariza- 
tion. SPOrt will observe a sky area with declination \S\ < 51.6° (~ 80% of 
the whole sky). The lifetime of the experiment is at least 1.5 years, but an 
extension is not unlikely; data will be binned in a number of pixels ranging 
from 662 to several thousands, depending on the effective duration of the flight 
and efficiency 

In this analysis we shall assume that pixels are distributed according to the 
HEALPix 5 package, with N stde = 8 or 16, and smoothed with a Gaussian 
beam of FWHM of 7°. The two choices of N stde allow to confirm that the 
best exploitation of the observational output is achieved when the average 
angular distance between pixel centers is ~ half FWHM. For the low (high) 
resolution cases we shall therefore have, on the whole sky, 768 (3072) pixels, 
whose centers lie at an average angular distance of ~ 7.3° (3.7°). Once polar 
caps are excluded, there remain 600 (2448) pixels, providing measures of the 

5 http : //www . eso . org/science/healpix/ 
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Stokes parameters Q and U. 



SPOrt polarimeters will provide no anisotropy data. Such information will be 
however available through other experiments (COBE, MAP) and we shall as- 
sume here that no peculiar problem arises in correlating SPOrt and anisotropy 
data. Of course, anisotropy data probe a different sky QXGcl. 1. e. the whole sky, 
excluding the area with declination \5\ < 20°, where galactic contamination 
is severe. (On the contrary, such contamination can be assumed to be under 
control, for polarization data; see, e.g., Bruscoli et al. 2002; Tucci et al. 2002.) 
For the coarse (fine) pixelization cases, as above, in this sky area there are 480 
(1984) pixels. For 376 (1360) of them, both anisotropy and polarization data 
are supposed to be available. 

Random noise (<Jp ix ) will be included in artificial data, assuming it to be 
uncorrelated, both among pixels and between Q and U, in the case of polar- 
ization measures. We shall assume = 1 /j,K (2//K) for temperature data (as 
expected for MAP measurements scaled to the SPOrt resolution). On the con- 
trary, one of the main variables we consider is the noise level for polarization 
data. 

SPOrt data will be characterized by a^ ix ~ 5.3 N ^ de \J^y^t /^K per pixel, 
where A is the experiment's lifetime in years while e indicates the detection 
efficiency Our analysis will consider cr^ x 's in the interval comprised from 0.5 
to 5 fxK per pixel, with coarse pixelization, corresponding to 1 to 10 fjK per 
pixel, with fine pixelization. 

The expected angular spectra Cf ,E,B ' x , were produced with the program 
CMBFAST 6 , starting from a spatially flat cosmological model with Vl m = 
0.35, Qb = 0.05, h = 0.65. Dark energy is an ordinary cosmological constant. 
Several values of n s were considered, but we shall report results only for n s — 1. 
Other close n s values yield quantitatively similar results and the potentiality 
of the experiment does not depend on n s , in any appreciable way. On the 
contrary, our analysis aimed to exploring r values ranging from 0.035 to 0.25. 
Cosmic variance was taken into account by considering 1000 independent re- 
alizations for each r and a^ ix pair. 



4 Data analysis 

The first aim of our analysis amounts to finding the likelihood distribution on 
the plane T-n s . The number of pixels for anisotropy and polarization (Nt and 
N P ) in our (artificial) data are different. Let then Tj be the anisotropy mea- 

6 http : //physics . nyu . edu/matiasz/CMBFAST/cmbf ast . html 
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sured in AT T pixels and Qj and Uj the Stokes parameters measured in Np pixels. 

In general, let us define vectors x = (T 1; , T Nt , Qi, , Qn p , U±, , Un p ), 

of N s = Nt + 2Np components, defining an observed state of anisotropy 
and polarization. Once a pair of values r-n s is assigned, the angular spectra 
Cf = (C; T , Cf , Cf) are univocally determined. On the contrary, a data vec- 
tor d, of iV s components, built from them, is just a realization of such model: 
once the N s component vector d is assigned, the pair of values r-n s is not 
univocally fixed. 

A function 

L(d\Cf) oc [detM]-5exp[- id T M-M] (1) 

shall then be built, to yield the likelihood of a given set of Cf 4 (i.e., of a pair 
of T-n s values), if d is observed. The main ingredient of L is the correlation 
matrix My = (xfxj) = Sjj + N^-; here Sy is the signal term and N^- is 
due to the noise. The components yield the correlation between the ith 
and jth elements of data vectors x corresponding to particular choices of Cf, 
i.e. of T-n s values. The construction of the (model dependent) signal term, 
however, does not require to build explicitly the vectors x. The procedure 
to be followed, in the case when both anisotropy and polarization data are 
available, is explicitly reported by Zaldarriaga (1998). The construction of 
the noise term is simpler, as we expect no noise correlation, and the matrix 
Nij = Sij o- T , P ix 2 (for i = 1, N T ) and = a PtPix 2 (for % = N T + 1, N s ) 
is diagonal. 

In what follows, the technical role of the Cf spectra will not need to be 
further outlined and the likelihood function will be explicitly considered to 
depend on r and n s . In particular, we assume that the most probable r-n s 
pair, for a given anisotropy-polarization state d, is the one which maximizes 
the likelihood. On the r-n s plane, we can then define curves connecting points 
where L corresponds to a given fraction of its top value. The area enclosed 
by each such curve corresponds to a given confidence level. In the figures 
below, such fractions were selected in order that the areas enclosed by the 
curves correspond to 68 %, 90 % and 99 % confidence levels. For the sake of 
simplicity, we shall refer to them as 1-2-3 a levels. 

Because of the role played by low-/ harmonics in our analysis, as above out- 
lined, cosmic variance may be quite significant. Henceforth, 1000 realizations 
of each model will be considered; 1-2-3 a curves will be obtained by averag- 
ing among the results of such realizations, so obtaining an average likelihood 
distribution. Of course, this is not sufficient to evaluate the expected result 
of a given experiment. We shall therefore implement average results with his- 
tograms, describing the distribution of results in various realizations, suitably 
defined for the various cases. 
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5 Results 



Because of the large number of realizations required, we considered likelihood 
distributions both for coarse and fine pixelization, corresponding to distances 
between pixel centers of ~ 7.3° or 3.7° (see section 3). A comparison of such 
distributions, for the cases we treated in both ways, confirms that, in general, 
a finer pixelization allows a better exploitation of data. Finding some general 
trends, however, does not require fine pixelization analysis. 

In Fig. 3, we report the likelihood distributions for r = 0.05 and 4 levels of 
a pix (f° r coarse pixelization). In Fig. 4, we report the likelihood distribution for 
the same r value and a^ ix = 2 /iK (for fine pixelization, corresponding to 1 /iK 
with coarse pixelization). Besides of the solid curves obtained by considering 
both anisotropy and polarization, dashed curves are shown, obtained from 
polarization data only. Stompor limits are explicitly shown in Fig. 4, while 
in Fig. 3 we report Stompor's 95 % confidence limit only. See the caption of 
Fig. 4 for further details on Stompor limits. 

From these figures it is clearly visible that upper limits on r can be significantly 
improved, in respect to previous analyses, for a^ ix < 2-3 /iK. This is true even 
if only polarization data are used. The main improvement obtained, when 
using anisotropy data as well, concerns n s , for which both lower and upper 
limits are obtained. When passing from coarse to fine pixelization, n s is even 
better constrained. If a^ ix is at the lowest levels considered, n s determination 
is better than in Stompor's case. 

For what concerns lower limits on r, the situation is more critical. Detecting 
a lower limit, however, should be the main aim of such kind of experiments. 
Fig. 4 shows that a lower limit, in average, can be detected at the 1-a con- 
fidence level at most. This is a case, however, when average estimates them- 
selves are insufficient. In Fig. 5, therefore, we show the distribution of 1-a 
lower limits found, with or without taking into account anisotropy, with fine 
pixelization. In the former case, we see that, in ~ 78 % of lower limit 

is found. Such percentage decreases to ~ 70 %, if no anisotropy data could be 
considered. For the sake of comparison, let us report that the above figures 
decrease to ~ 76 % and ~ 64 %, with coarse pixelization. 

Going below a noise level of 2 /iK per pixel (with fine pixelization), although 
highly problematic from an experimental point of view, would be however 
decisive. This is shown by Fig. 6, where we see that, for a^ ix = 1 /iK, even 
the average 3-cr curve indicates a detection of t > 0. If the distribution of 
lower limits in various realizations is considered, we find that, only in ~ 2 % 
of cases, no 3-cr lower limits would be detectable. 

Let us now consider the cases with greater r values. They are reported in 
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Fig. 7 and 8, for r = 0.10 and 0.15. These figures are all obtained with coarse 
pixelization. They confirm that significant improvements of Stompor upper 
limits on r are obtainable for a^ ix < 2-3 /iK. Moreover, lower limits become 
safer than in the r = 0.05 case. Even with a^ ix = 2 /jK, if r ~ 0.15, a \-a 
lower limit is likely to be obtained and a 2-cr lower limit might be obtained 
as well. This is an important point, as such noise levels are not far from those 
obtainable, in the best possible conditions, by the SPOrt experiment. 

It is also significant to go above average values, for lower limits, in all previous 
cases. In Fig. 9 to 11, we report the fractions of realizations, for which lower 
limits were detected, both using anisotropy-polarization data and polarization 
only. Among various possible comments, these plots confirm that, in order to 
detect r (unless > 0.15), the noise level of the SPOrt experiment is too high. 
It may also be significant to compare these figures with the expected noise 
level in the PLANCK measurements of polarization. Weighting the various 
channels in order to minimize variance, it is expected that a^ ix ~ 9 fiK, for a 
FWHM of 7' (see, e.g., Hu & Okamoto 2002). If a scaling a^ ix oc angle^ 1 is 
assumed, this yields a^ ix ~ 0.4 yuK, if reported to our fine pixelization. At such 
noise level, no magnification of the intrinsic cosmic variance seems to occur. 



6 Conclusions 

In this paper we performed a likelihood analysis, to determine how a polariza- 
tion experiment could break the degeneracy between n s and r determinations, 
found in pure anisotropy analyses. To this aim we built artificial data, consis- 
tent with an experiment in progress to detect polarization, but pushing noise 
even well below the expected experimental level. The experiment considered 
will cover 80 % of the whole sky and artificial data were built assuming such 
sky coverage. If data were available for the full sky, a slight improvement of 
the signal can be expected. Artificial data were then analyzed. In respect to 
future observers, we have however the critical advantage to be able to analyze 
as many skies (realizations) we need. This is significant, in our case, as the 
n s -r degeneracy, in principle, is readily overcome from low-/ spectral compo- 
nents, where, however, cosmic variance can be critical. Our results, however, 
indicate that cosmic variance is not such a severe limitation, if a sufficiently 
low noise level is attained. This is one of the most important conclusions of 
our analysis, as it confirms that, to detect r, polarization experiments should 
aim to high sensitivity, while high resolution does not bear a direct relevance. 

The critical issue, in this respect, is whether and when firm lower limits on 
r are obtainable. The noise level to be attained, to implement such aim, ob- 
viously depends on the physical value of r itself. Available constraints on r, 
besides that from CMB anisotropy data, can be also related to the detection of 
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hydrogen lines in high-z QSO spectra. This however leaves us a still significant 
range of possible cosmic opacities to CMB photons. 

If r is not negligibly low, the noise level below which we expect that lower lim- 
its on r can be obtained is a^ ix ~ 1-2 /xK, for an instrument with an angular 
resolution similar to COBE. A detection, at higher a^ x , cannot be however 
excluded. We compared such general conclusion with the actual expected fea- 
tures of the SPOrt experiment and found that, for some of the greatest values 
for r allowed by the cosmic reionization physics, and assuming the best pos- 
sible performance of the experiment, a detection of r, at the 2-a level, is 
possible. 
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Fig. 2. Dependence of r on the reionization redshift z r i in various cosmological 
models. The only model with total density parameter Otot 7^ 1 is the second one. 
In the third model, dark energy is due to a scalar field, self-interacting through a 
Ratra-Peebles potential (1988), with A RP /GeV = 10 3 . In all models n b h 2 = 0.022. 
We assume (almost) all baryons to be ionized after z r %. 
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Fig. 3. Average likelihood distributions, in a model with r = 0.05, if the pixel-pixel 
angular distance is 7.3°. Solid curves are obtained considering both anisotropy and 
polarization data; dashed curves are obtained with polarization data only; the cor- 
responding confidence levels are reported in the text. 
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Fig. 4. Average likelihood distributions, in a model with r = 0.05 and o^ ix = 2 fiK, 
if the pixel-pixel angular distance is 3.7°. Solid and dashed curves are as in Fig. 3. 
The limits found by Stompor at al (2001), using anisotropy data from balloon 
experiments, are also shown: The dot-dashed contour limits the 95 % confidence 
area, without any prior (this contour is reported in other plots as well). Assuming 
that n s < 1.2, they set the limit r < 0.4, at the 95% c.L, (long-short dashed 
contour). With the prior n s = 1, instead, r < 0.26, as shown by the arrow followed 
by a dotted line. Without any prior, however, they find that, at the 99 % c.L, n s > 0.8 
(long dashed line). 
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Fig. 5. Distribution of 1— a lower limits, for r = 0.05 and o^ ix = 2 fiK, if the 
pixel-pixel angular distance is 3.7°. The l.h.s. (r.h.s.) plot includes (excludes) 
anisotropy data. 
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Fig. 6. Average likelihood distributions, in a model with r = 0.05 and a^ ix = 1 fiK, 
if the pixel-pixel angular distance is 3.7°. Solid and dashed curves are as in Fig. 3. 
This plot shows that, at this noise level, r could be fairly reliably detected. 
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Fig. 7. Average likelihood distributions, in a model with r = 0.10, if the pixel-pixel 
angular distance is 7.3°. Solid and dashed curves are as in Fig. 3. 
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Fig. 8. Average likelihood distributions, in a model with r = 0.15, if the pixel-pixel 
angular distance is 7.3°. Solid and dashed curves are as in Fig. 3. 
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Fig. 9. Fractions of realizations for which lower limits on r are detected, at the 
l-o" confidence level, for different r and a^ ix values. Solid lines are obtained with 
anisotropy and polarization data. Dashed lines are obtained just with polarization. 
The values of o^ ix reported in the plot, refer to a coarse pixelization analysis. For 
low values of n/Nt Q t, random noise causes some irregular behaviour. 
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